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Abstract 

This paper is concerned with the study of insurance related derivatives on financial mar- 
kets that are based on non-tradable underlyings, but are correlated with tradable assets. We 
calculate exponential utility-based indifference prices, and corresponding derivative hedges. 
We use the fact that they can be represented in terms of solutions of forward-backward 
stochastic differential equations (FBSDE) with quadratic growth generators. We derive the 
Markov property of such FBSDE and generalize results on the differentiability relative to 
the initial value of their forward components. In this case the optimal hedge can be repre- 
sented by the price gradient multiplied with the correlation coefficient. This way we obtain 
a generalization of the classical 'delta hedge' in complete markets. 
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Introduction 

In recent years more and more financial instruments have been created which are not derived from 
exchange traded securities. For instance in 1999 the Chicago Mercantile Exchange introduced 
weather futures contracts, the payoffs of which are based on average temperatures at specified 
locations. Another example of derivatives with non-tradable underlyings are catastrophe futures 
based on an insurance loss index regulated by an independent agency or simply derivatives based 
on equity indices such as S&P or DAX. 

Financial or insurance derivatives of this type are impossible to perfectly hedge, since it is 
impossible to trade the underlying variable that carries independent uncertainty. To circumvent 
this problem, in practice one looks for a tradable asset that is correlated to the non-tradable 
underlying of the derivative. Even though investing in the correlated asset cannot provide a total 
hedge of the derivative, and a non-hedgeable basis risk remains, it is better than not hedging at 
all. 

In the following we will investigate utility-based pricing principles for derivatives based on 
non-tradable underlyings. Moreover we will show how the derivatives can be partially hedged by 
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investing in correlated assets. We present explicit hedging strategies that optimize the expected 
utility of a portfolio of such derivatives. To this end we will establish some structure and 
smoothness properties of indifference prices such as the Markov property and differentiability 
with respect to the underlyings. Once these properties are established, we can explicitly describe 
the optimal hedging strategies in terms of the price gradient and correlation coefficients. This 
way we obtain a generalization of the classical delta hedge of the Black-Scholes model. 

The hedging of claims based on non-tradable underlyings has already been studied by many 
authors, see for example |HH02| . [Hen02j . |MZ04| . |Dav06| . |Mon04j . |AIP07| . As a common 
feature of all these papers, optimal hedging strategies are derived with standard stochastic 
control techniques. The essential components of this analytical approach consist in a formulation 
of the optimization problem in terms of HJB partial differential equations, and the use of a 
verification theorem and uniqueness result in order to obtain a representation of the indifference 
price and the optimal control strategy. We instead employ an approach with a stochastic focus. 
It starts with the well-known observation that the maximal expected exponential utility may be 
computed by appealing to the martingale optimality principle which leads to a description of 
price and optimal hedging strategy in terms of a forward-backward stochastic differential equation 
(FBSDE) with a nonlinearity of quadratic type (see [REKOOj . |HIM05]). This immediately 
implies that the utility indifference price resp. hedge is equal to the difference of initial states 
resp. control processes of two FBSDE with a quadratic nonlinearity in the generator. The 
forward component is given by a Markov process describing the non-tradable underlying. The 
main mathematical contribution of this paper is that it provides simple sufficient conditions for 
general FBSDE with quadratic nonlinearity to satisfy a Markov property, and - for the BSDE 
component - to be differentiable with respect to the initial condition of the forward equation. 
The techniques for proving differentiability of BSDE with quadratic nonlinearity have been 
developed independently in [BC07| and [AIDR07] . Unfortunately, the setup of both papers 
is not general enough to cover the BSDE needed to calculate exponential indifference prices. 
Therefore, a slight generalization of these differentiability results is given in the last section of 
this paper. 

As a consequence of the explicit description of indifference prices and hedges in terms of 
the solution processes of the FBSDE, and in view of the smoothness results mentioned, it is 
straightforward to describe optimal hedging strategies in terms of the indifference price gradi- 
ent and the correlation coefficients explicitly. An economics related contribution of the paper 
is that the framework presented allows to refine the results obtained for example in [MZ04] . 
|Dav06j . Firstly, no longer we need to impose any restrictions on the coefficients of the diffusion 
modeling the tradable asset price. More importantly, the BSDE techniques allow to deal with 
multidimensional underlyings and traded assets. In the approach based on the HJB equation, 
a solution of the PDE is obtained by using an exponential Hopf-Cole transformation that in 
general seems to require that there exists only one traded asset. In practice many derivatives 
are based on more than one underlying, such as spread options or basket options. In order to 
illustrate how to hedge with more than one asset, we will study in more detail so-called crack 
spreads, which are written for instance on the difference of crude oil futures and kerosene prices 
(see Example [L2] and 0~§} . 

Finally we address the pricing of derivatives by the marginal utility approach. If a company 
wishes to trade risk not covered by securities on an exchange, they are forced to go outside the 
exchange to get tailored products to serve their specific needs. These deals that do not go through 
the exchange trading (although the underlyings may be traded there) and are done directly 
between buyer and seller are called over-the-counter (OTC). For example, airlines regularly 
make this kind of OTC deals in order to protect themselves against kerosene price fluctuations, 
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which underlines that the amount of money involved in this type of deals is non-negligible! 
Investment banks offering OTC deals face the problem of finding a fair price of these agreements. 
Indifference prices are often a reasonable solution. However, they are not linear! The standard 
way out, as suggested in the economics literature, is pricing by marginal utility. The marginal 
utility price is the differential quotient of the indifference price with respect to a marginal amount 
of the derivative. Here again the first thing to verify is the differentiability of the FBSDE. This 
in turn allows to derive the dynamics of the marginal utility price as a BSDE with a driver 
satisfying a random Lipschitz condition. 

BSDE with generators of quadratic nonlinearity in the control variable (which will in the 
sequel sometimes simply be called quadratic BSDE) are described by equations of the type 

Y t = £ + J T f(s, Y s , Z s )ds - J Z s dW s , < t < T, 

where / is a predictable function satisfying \ f(t,y,z)\ < C(l + \y\ + \z\ 2 ) with some constant 
C. Our differentiability results are based on the assumption that the derivative to be hedged, 
denoted by £, is essentially bounded. This guarantees that the integral process J Q ZdW is a so- 
called BMO martingale, and hence the density process of a new equivalent probability measure, 
say Q. By switching to the measure Q one can derive moment estimates needed in order to 
prove differentiability. The assumption that the derivative has to be bounded seems to be a 
disadvantage of using BSDE in the stochastic approach instead of working with the HJB partial 
differential equation in the analytical approach. In practice, this is of no importance. 

The paper is organized as follows: in Section [T] we introduce the model, in Section [2] we 
briefly recall results from [HIM05] concerning the solution of the problem of exponential expected 
utility maximization in terms of stochastic control problems and FBSDE with nonlinearities of 
quadratic type. In Section [3] we show structure properties of indifference prices of derivatives 
based on a non-tradable Markovian index process. In Section 2] we derive explicit formulas 
for the optimal hedges of such derivatives, and in Section [5] we describe the dynamics of the 
marginal utility price. All the economics related results are based on mathematical properties 
of quadratic FBSDE, which will be proved in the last section. 



1 The model 

Let d S N and let W be a (i-dimensional Brownian motion on a probability space (fi, T, P). We 
denote by (Ft) the completion of the filtration generated by W. Suppose that a derivative with 
maturity T > is based on a R m -dimensional non-tradable index (think of a stock, temperature 
or loss index) with dynamics 

dR t = b(t, R t )dt + p(t, R t )dW t , (1) 

where b : [0, T] x R m — > R m and p : [0, T] x R m -> R mxd are measurable deterministic functions. 
Throughout we assume that there exists a C € M+ such that for all t & [0, T] and x, x' E R m 

\b{t,x)-b(t,x')\ + \p{t,x)-p{t,x')\ < C\x-x'\, 
[ ' \b{t,x)\ + \p{t,x)\ < C{\ + \x\). 

We consider a derivative of the form F(Rt), where F : R m — > R is a bounded and measurable 
function. Note that at time t, the expected payoff of F(Rt), conditioned on Rt = r, is given by 
F{R t rf), where R l ' r is the solution of the SDE 

Rtf = r + T b(u, B*> r )du + f p(u, Rf)dW u , s G [t, T]. (2) 
Jt Jt 
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Our correlated financial market consists of k risky assets and one non-risky asset. We use 
the non-risky asset as numeraire and suppose that the prices of the risky assets in units of the 
numeraire evolve according to the SDE 

dS\ = Si(an(t, R t )dt + (3i(t, R t )dW t ), i = l,...,k, 

where oti(t, r) is the ith component of a measurable and vector- valued map a : [0, T] x R m — » R fc 
and @i(t,r) is the zth row of a measurable and matrix- valued map : [0, T] x W 11 — ► R fcxd . 
Notice that W is the same R d -dimensional Brownian motion as the one driving the index process 
(flj). and hence the correlation between the index and the tradable assets is determined by the 
matrices p and (3. 

In order to exclude arbitrage opportunities in the financial market we assume d > k. For 
technical reasons we suppose that 

(Ml) a is bounded, 

(M2) there exist constants < e < K such that el k < (P(t,r)fl*(t,r)) < KI k for all (t,r) € 
[0, T] x R m , 

where /3*(t,r) is the transpose of f3(t,r), and I k is the ^-dimensional unit matrix. 

Before we proceed with the model description we will illustrate the range of possible appli- 
cations by giving some examples of derivatives our model may apply to. 

Example 1.1. Weather derivatives are typical example of financial instruments derived from 
non-tradable under lyings. One of the most common types of weather derivatives are based on 
so-called accumulated heating degree days (cHDD). The heating degree of a day with average 
temperature r in Celsius degrees is defined as HDD = max{0, 18 — r}, i.e. HDD describes 
the (positive) difference between the average daily temperature measured and the temperature 
above usually rooms are heated. The cHDDs are defined as a moving average sum of HDDs over 
a fixed time length, for instance a month. Real data shows the cHDD to be almost lognormally 
distributed, and therefore they can be modelled as geometric Brownian motions (see [DavOlj ). 
This means that in ([T]) we would have to choose b(t, Rt) = a±Rt and p(t, Rt) = a^Rt, with a\ G R 
and OLi £ R \ {0} depending on the season. Tradable assets that are more or less correlated with 
average temperatures are for example electricity futures and natural gas futures. 

The derivative explained in the next example is based on more than one underlying. 

Example 1.2. Spread options in general involve two or more underlying structures (prices, 
indices, interest rates and many other possible quantities), and measure the distance between 
them. We do not go into details since spread options are well-known (see |CD03| for an overview). 
For simplicity we refer to a 2 dimensional example of Crack spreads. 

Crack spreads consist in the simultaneous purchase or sale of crude against the sale or 
purchase of refined petroleum products. We concentrate on the kerosene crack spread, which 
pits crude oil price (co) against kerosene price (ke). A company producing kerosene (from crude 
oil) wishes to cover part of its risk arising from a sudden boost of the crude oil price by buying 
kerosene crack spreads. It thereby faces the problem that kerosene trading is not done on a 
sufficiently liquid market to warrant a futures contract or some other type of exchange-traded 
contract. So derivative contracts of this type must be arranged on over-the-counter basis. 

Knowing that the price of heating oil (ho) is highly correlated with the kerosene price - except 
during the Iraq war - crack spreads themselves can be hedged by using heating oil futures. 
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We model prices in the following way, where the superscripts represent the underlying prod- 
ucts, 

dR k t e = R k t e (hdt + l2 dWl + jidWf + -f 4 dWf) 

dR c t ° = R c t ° (b 2 dt + jidWl) 

dS?° = S?° (b 3 dt + fcdW? + /3 2 dW?) , 

where we assume that bi,b 2 ,bs S 71)72)73) 74) fhifh € K\{0} and the correlation between 
heating oil and kerosene is given by a = (72/^1 + ^3^2) /\/ (t| + 73 + + P 2 )- 

A European call on the spread is of the form £,(R^ e , Sf) = (R^ e - - K) + , with K being 
the strike. 

Throughout let U be the exponential utility function with risk aversion coefficient ij > 0, i.e. 

U{x) = -e~^ x . 

In what follows let (t, r) G [0, T] x M. m . By an investment strategy we mean any predictable 
process A = (A*)i<j<fc with values in M fc such that the integral process J * \ l r ^f is defined for all 

i € {1, . . . , fc}. We interpret A* as the value of the portfolio fraction invested in the i-th asset. 
Investing according to a strategy A leads to a total gain due to trading during the time interval 
[t, s] which amounts to G^' 1 = Y^=\ It ^u^§~- We will denote by G^ ,t,r the gain conditional on 
Rt = r. 

Remark 1.3. As one can see the wealth process is given by 

G s' t,r = E f X ul a *( u > Rt u)du + (3i(u, R l /)dW u )], 
i=i J* 

and hence does not depend on the value of the correlated price process! This is a feature of our 
model that will later imply the indifference price at time t to depend only on the value of the 
index process at a given time t. 

Let A 1 ' 1 " be the set of all strategies A such that E \X s /3(s, R t s ' r )\ 2 ds < 00 and the family 

{e~ vGr : t is a stopping time with values in [t, T]} is uniformly integrable. If A G A t,r , then 
we say that A is admissible. We use the same admissibility criteria as in Section 2 in |HIM05j . 
so that later we may invoke their results. The maximal expected utility at time T, conditioned 
on the wealth to be v at time t and the index to satisfy Rt = r, is defined by 

V°(t, v, r) = su V {EU(v + G^' r ) : A G A l ' r }. (3) 

One can show that there exists a strategy n, called optimal strategy, such that EU(v + G^,' i,r ) = 
V°(v,t,r). The convexity of the utility functions implies that -zr is a.s. unique on [t,T], and it 
follows from Theorem 7 in |HIM05| that ir G A 1 ' 1 '. 

Suppose an investor is endowed with a derivative F{Rt) and is keeping it in his portfolio 
until maturity T. Then his maximal expected utility is given by 

V F (t, v, r) = su V {EU{v + G^' r + F(R^f)) : A G A*'"'}. (4) 

Also in this case there exists an optimal strategy, denoted by 7?, that satisfies EU(v + G 1 ^ *' r + 
F{R^ r )) = V F {v,t,r). 
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The presence of the derivative F(Rt) leads to a change in the optimal strategy from tt to tt. 
The difference 

A = TT — TT 

is needed in order to hedge, at least partially, the risk associated with the derivative in the 
portfolio. We therefore call A derivative hedge. In the following sections we shall analyze by 
how much the optimal strategies change if a derivative is added to the portfolio, and we aim at 
getting an explicit expression for the derivative hedge A. 

One can easily show that for all (£, r) G [0, T] x R m there exists a real number p(t, r) such 
that for all v G R 

V F (t,v-p(t,r),r) = V°(t,v,r). 

If an investor has to pay p(t, r) for the derivative F(R^), then he is indifferent between buying 
and not buying the derivative. Therefore the number p(t, r) is called indifference price at time 
t and level r. 

It turns out that the derivative hedge A is closely related to the indifference price of the 
derivative. The derivative either diversifies or amplifies the risk exposure of the portfolio. The 
difference between tt and tt measures the diversifying impact of F. The price sensitivity, i.e. the 
derivative of p relative to the index evolution, is also a measure of the diversification of F (which 
will be called diversification pressure of the derivative F). We will see that the derivative hedge 
is indeed equal to the price sensitivity multiplied with some correlation parameters. 

The problem of finding the optimal strategies tt and tt is a standard stochastic control 
problem. One can tackle it by solving the related HJB equation, using a verification theorem 
and proving a uniqueness result. This approach has been chosen for example in |AIP07j . Here, 
however, we prefer a stochastic approach, using the fact that the stochastic control problem can 
be solved by finding the solution of a backward stochastic differential equation (BSDE). In the 
following section we briefly recall the definition of a BSDE. 



2 Solving stochastic optimal control problems via BSDE 

Let 7i 2 (M. d ) be the set of all Revalued predictable processes £ such that E \(t\ 2 dt < oo, and 
let 5 2 (R) be the set of all R- valued predictable processes 5 satisfying E ^sup se [ 0) T] |£s| 2 ) < °°- 
By <S°°(R) we denote the set of all essentially bounded R- valued predictable processes. Let £ 
be .Fr-measurable and / a predictable mapping defined on 17 x [0, T] x R x R d with values in 
R. A solution of the BSDE with terminal condition £ and generator f is defined to be a pair of 
processes (Y, Z) G S 2 (R) x H 2 (M. d ) satisfying 



Y t = £- f T Z s dW s + f T f(s,Y s ,Z s )ds. 
Jt Jt 



Let us now come back to our control problem of finding the optimal investment strategy tt and tt 
respectively. It is known that there exists a quadratic BSDE which solves these control problems 
(see for example [HIM05J). We first specify the generator of the suitable BSDE, starting with 
tt. 

Fix again (t,r) G [0,7] x R m . Let 0(t,r) = P*{t,r)(p(t,r)P*(t,r))- l a(t,r) and C(t,r) = 
{x(3(t,r) : x G R fe }. Observe that our assumptions imply that $(t,r) is bounded. The distance 
of a vector z G R d to the closed and convex set C(t,r) will be defined as dist(z, C(t, r)) = 
minjjz — u\ : u G C(t,r)}. Let / be the deterministic function 

/ : [0,71 xl m xl^K, (t,r,z) ^ zti(t,r) + —\m,r)\ 2 - ^dist 2 (z + r), C(t, r)). 

2rj 2 r\ 
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Since d > k, we have to find the orthogonal projection of the d— dimensional vector z to the 
linear space C(t,r) of image strategies. In |HIM05j the set C(t,r) is understood as imposing 
restrictions on the investor when trading in the market that happen to be convex in the setting 
given. 

Notice that / is differentiable in z and satisfies the growth condition 

\f(t,r,z)\<c(l + \z\ 2 ) a.s. 

with some c G M+. The growth condition guarantees that there exists a unique solution 

(y*' r , Z l ' r ) G «S°°(M) &> H 2 (R d ) of the BSDE 

= F(R% r ) - f Z\fdW u - f f(u, Z^du, s G [t, T], (5) 

J s J s 

(see Theorem 2.3 and 2.6. in [KobOOj). Notice that the terminal condition of the BSDE stems 
from a standard forward SDE. The system of equations consisting of ([2]) and ([5]) is often called 
forward-backward stochastic differential equation (FBSDE). 

The conditional maximal expected wealth, or in other words the value function of our stochas- 
tic control problem, is equal to the utility of the starting point of the BSDE, i.e. 

V F (t,v,r) = -e-^-^'") 

(sec Theorem 7 in [HIM05 ). Moreover we can reconstruct the optimal strategy tt starting from 
Z. To this end denote by T\-c(t,r){ z ) t ne projection of a vector z € M d onto the linear subspace 
C(t,r). If Rt = r, then the optimal strategy ■nt on [t, T] satisfies 

t? s (}{s,RY) = n c ^ RY) [ZY + U{s,RY% s G [t,T]. (6) 

The last statement follows equally from Theorem 7 in [HIM05] . 
Analogously, let (Y t,T ', Z t,r ) be the solution of 

Ys' r = ~ ^ Z' u r dW u - f f(u, R 1 /, Z^du, s G [i, T], (7) 

J s J s 

which represents a stochastic control problem as above, just without the derivative as terminal 
condition i.e. the derivative is not in the portfolio. In this case the maximal expected utility 
verifies 

V°(t,v,r) = - e -^ v - y i' r \ 
and the optimal strategy 7r on [t, T] satisfies 

TTsP&Rtf) = n c ^ RY) [ZY + U{s,RY% s G [t,T]. (8) 

Since 11^^ R t, r ^ is a linear operator, the derivative hedge is given by the explicit formula 

A s p(s,Rln = n c[s ^ r) [Z^ - Z'/}, 
which will be further determined in the subsequent sections. 
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3 The Markov property of the indifference prices 



In this section we will establish the Markov property of the indifference prices. This will follow 
from the fact that the solutions of the BSDEs © and (|7j) are deterministic functions of time and 
the underlying. To give the precise statement we need to introduce the following cr-algebras. Fix- 
ing t G [0, T], we denote by T> m the cr-algebra generated by the functions r i— > E[J^ <f>(s, R t s ' r )ds], 
where t G [0, T] and is a bounded continuous R— valued function. 

Moreover we assume that the mapping (t, r) i— > $(t, r) be Lipschitz continuous in r, noting 
that due to (Ml) and (M2) this is guaranteed if (3 and a are Lipschitz continuous. 

Lemma 3.1. There exist £>[0,T] ® T> m -measurable deterministic functions u and u : [0,T] x 
R m -> R such that 

Y l s ' r = u(s,Rl' r ) and Yg' r = u(s, R 1 /'), 
for P <g> X-a.a. (uj, s) G f2 x [t, T] . 

Proof. The generator function / is a polynomial of the components of z of at most second degree. 
This implies, together with the assumption that i? is Lipschitz continuous in r, that there exist 
functions f n : [0, T] x R m x M. d — > R, globally Lipschitz continuous in z, such that for all compact 
sets K C R m x R rf the sequence f n converges to / uniformly on [0, T] x K. Thus the statement 
follows from Theorem 16.61 □ 

Lemma 13.11 immediately implies that there exists a nice version of the indifference price p as 
a function of (t, r). 

Theorem 3.2. There exists a B[0,T]®T> rn -measurable deterministic function p : [0,T]x]R m — > R 
such that for all v G R, (t, r) G [0, T] x R m 

V F {t,v -p(t,r),r) = V°(t,v,r). (9) 

Proof. Let v G R, (i, r) G [0,T] x M m be given. Recall that 1/ F (u,t,r) = -e-''( 1 '- y t t ' r ) and 
V°(v,t,r) = —e~ r] ^"~ Yt Then put p(t,r) = u(t,r) — u(t,r), where u and u are given from 
Lemma 13. 11 □ 

In the remainder the function p is always assumed to be measurable in both t and r. In fact 
it inherits this property from the functions u and u. 

We now turn to an explicit description of the optimal strategies, and in particular their 
difference, the derivative hedge. These will be derived from the BSDE solutions of the preceding 
section. We start by noting that similarly to the indifference price the optimal strategies only 
depend on the time and the index process R. 

Theorem 3.3. There exist B[0,T] <8> T> m -measurable deterministic functions v and v, defined 
on [0, T] x R m and taking values in R d such that for (t,r) G [0, T] x R m , the optimal strategies, 
conditioned on Rt = r, are given by tt s = is(s, R^ r ) and tt s = u(s, R l s T ) for all s G [t, T]. 

Proof. Fix (t, r) G [0, T] x R m . Theorem 16.61 implies that there exist B[0, T] ® X> m -measurable 
deterministic functions v and v mapping [0, T] x R m to R m such that for all s G [t, T] 

Z l f = v(8,R?)p{s,I% r ) and ^ = 5( g ,^)p( g ,l§ r ). 

Now let 7 (t,r) = U c{tir) [v(t,r)p(t,r) + ^(t,r)} and j(t,r) = U c{tir) [v(t,r)p(t,r) + ±0(t,r)]. 
Then, by © and ([8]), the optimal strategies conditioned on Rt = r satisfy 

9 s f3(s, RY) = 7( a , RY) and 7r s f3(s, 2# r ) = 7 (a, 
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for alls £ [t,T]. Since the rank of (3(t, r) is k, then both V(t, r) = ^(t,r)f3*(t,r)((3(t,r)(3*(t,r))~ 1 
and u(t, r) = j(t, r)(3*(t, r)(/3(i, r)f3*(t, r))" 1 are well defined. Then uniqueness of it and tt yields 
the result. □ 

Remark 3.4. Theorem 13 . 31 implies that the optimal strategies are the so-called Markov controls. 

We close this section by noting that Theorem 13.21 implies a dynamic principle for the indif- 
ference price. Abbreviate A = A°' r for some r £ W 71 . For any stopping time r < T and J- T - 
measurable random variable G T let V (r, G r ) = esssup {E[U(G T + G$ T + F(B% r )\T T ] : A e A}. 
Similarly we define V°(r, G T ). 

Corollary 3.5. We Ziawe 

V F (T,G T -p(T,R° T n) = V°(T,G T ). 
Proof. As is shown in Prop. 9 in [HIM05J, the value function V F satisfies the dynamic principle 
V f (t,G t - P (r,R° T n) = U(G T -p(T,R /) - Y T °> r ). 

Since p{T,R° T ' r ) = Y?'^ - Y^ = Y?' r - Y r °' r we obtain V f (t,G t -p(r,i# r )) = U(G T - 
Y T °> r ) = V°(T,G T ). □ 

4 DifFerentiable indifference prices and explicit hedging strate- 
gies 

If we impose stronger conditions on the coefficients of the index process R and the function F, 
then we can show that the price function p is differentiable in r, and we can obtain an explicit 
representation of the derivative hedge in terms of the price gradient. To this end we need to 
introduce the following class of functions. 

Definition 4.1. Let n, p > 1. We denote by B nxp the set of all functions h : [0, T] x R m — > 

M nxp , (t, x) i — > h(t,x), differentiable in x, for which there exists a constant C > such that 



sup (M; ) e [ ,T]xR™ YT=i nip - c > for a11 1 e [°' T ] we have su p^6r™ J i+ff < c> and x 

9h Q^ is Lipschitz continuous with Lipschitz constant C. 

We will assume that the coefficients of the index diffusion satisfy in addition to (Rl) 
(R2) p £ B mxd , b £ B mxl , and 

(R3) F is a bounded and twice differentiable function such that 

VF-pe B 1 *" and J>(*,r)^F(r) + ~ 2 E bA^ft^M € ^ " 
«=i «,i=i * ^ 

The next result guarantees Lipschitz continuity and differentiability of the functions u and S 
obtained from Theorem 13. 11 

Theorem 4.2. Suppose that (Rl), (R2) and (R3) are satisfied. Besides, suppose that the 
volatility matrix f3 and the drift density a are bounded, Lipschitz continuous in r, differentiable 
in r and that for all 1 < i < k, 1 < j < d the derivatives V r /3jj and V r aj are also Lipschitz 
continuous in r. Then the functions u and u are Lipschitz continuous in r, and continuously 
differentiable in r. 



\h(t,x)\ 
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Proof. The theorem follows from Lemma [6. 3 1 and Theorem 16.71 in Section [6l All we have to show 
at this stage is that the assumptions of both results are satisfied. We only show Conditions (|20p 
and (|23p . since the remaining ones are easily seen to be fulfilled. 

Notice that the conditions on a and /3 imply that # is differentiable in r, and that i? and V r i? 
are globally Lipschitz continuous in r. Moreover, since i? is bounded, t? 2 is Lipschitz continuous 
in r, too. Recalling the definition of the generator /, note further that 

|dist 2 (z + -m, r), C(t, r)) - dist 2 (z + -0(t, r'), C(t, r))\ 

< 2(-||0|| oo + |z|) |dist(z + -ti(t,r),C(t,r)) - dist(s + -0(t, r'), C(t, r))| 
77 77 77 



< 2(i||tf||oc + M 
77 



I* + r) - Hc M (z + -^(t, r))| -\z + -&(t, r') - U c[Lr) (z + r'))| 



< 4(-\\&\\oo + \z\) - \#(t,r) -fi(t,r') I , te [0,r],r,/eR m ,z€R rf . 
77 77 

This shows that there exists a constant X € M+ such that r, z) — f(t, r' , z)\ < K(l + \z\)\r — 
r'\, and hence Assumption (I20p of Lemma 16.31 is satisfied. 

Observe that U c ^ t ^(y) = r) for all y G and hence the mapping r 1— > 

n C ( (r )(z + r)) is differentiable. Consequently, also / is differentiable and for t G [0, T], 
r£R m and 2 G R d we have 

V r /(t, r, z) = zVMt, r) + -&(t, r)V r 0(t, r) 

V 

-r, (z + ~ti(t, r) - U c(ttr) (z + ^{t, r))J Q V r 0(t, r) - V r n c(t>r) (z + r)) 

By using that i?, V r ?9, /3 and V r/ 9 are Lipschitz continuous and bounded, it is straightforward 
to show that for all t G [0, T], r, r' G R m and z, z' G R d we have 

I V r /(t, r, z) - V r /(t, r', z')\ < K(l + |z| + \z'\)(\r - r'\ + \z - z'\), 

and hence the generator satisfies Assumption (|23p of Theorem 16.71 

Now Lemma 16.31 yields the Lipschitz continuity in r of the functions u and u. Theorem 16.71 
implies the differentiability of Y t,r and Y t,r with respect to r, and hence also of u and u. □ 

As an immediate consequence we obtain smoothness of the indifference price function. 

Corollary 4.3. Suppose that the assumptions of Theorem \^.S\ are satisfied. Then the indifference 
price function p is continuously differentiable in r. 

Having shown smoothness of the indifference price, we can finally derive an explicit formula 
for the derivative hedge in terms of the price gradient. To this end we denote the conditional 
derivative hedge by A(t, r) = v(t, r) - u(t, r), (t, r) G [0, T] x R m . 



Theorem 4.4. Under the assumptions of Theorem \4-ty and with the notation of Section® the 
derivative hedge satisfies 

A(t,r) = -V r p(t,r)p(t,r)p*(t,r)(p(t,r)(3*(t,r))-\ (t,r) G [0, T] x R m . (10) 

Remark 4.5. Note that Theorem 14.41 implies that the derivative hedge at time t depends only 
on R t . 
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Proof of Theorem\4^ Note that C(t,r) is a linear subspace of R d for all (t,r) G [0,T] x M m . 
Therefore, the projection operator IIc(t,r) i s linear and hence 

A(i,r) = (T^O-T^r))/?*^^^^)/?*^^)- 1 

n c(iir) [Zf + jU(i,r)] - nc (tjr) [Z*' p + ^(t,r)A (?(t,r)<J){t,r)ir{t,r))- x 

It follows from Theorem 16. 71 that Z*' r — Z l t ' r = (V r u(i, r) — V r u(t, r)) p(t, r) = —V r p(t, r)p(t, r), 
and hence we obtain the result. □ 

If the market consists of only one risky asset, then the optimal strategy simplifies to the 
following formula. 

Corollary 4.6. Let k = 1. Then the derivative hedge is given by 

(P(t,r),V r p(t,r)p(t,r)) gtl ££i ^p(t, r)p 3 -,(t, r) 
( ' )_ " E^/fM ' (t ' r)Gi °' TJxM 

Proof. Fix (i,r) G [0, T] x R m . Note that C(t,r) = {x(3(t,r) : x G M} is a one-dimensional 

subspace of M d For all z = (z 4 )i<*<d e R d let 5 (z) = ^jg^ = ^^f^ . Then 

g(z)j3(t, r) is the orthogonal projection of z onto C(i, r). Thus Theorem 14.41 yields that A(t, r) = 
- 9 (V r p(t,r)p(t,r)). □ 

Remark 4.7. 

1) Suppose the derivative F{Rt) is traded on an exchange. By pretending the price observed 
is approximately equal to an indifference price, the hedging formula (jlOl) provides a very simple 
tool for hedging the derivative. Notice that the risk aversion coefficient ry does not appear 
explicitly in (jlOD . 

2) If k = d and the matrices j3(t, r) are all invertible, then our financial market is complete and 
the derivative F(Rt) can be fully replicated. Moreover the derivative hedge satisfies 

A{t,r) = -V r p{t,r)p{t,r)j3~ 1 {t,r). 

If S is chosen to be the index, i.e. R = S, then we obtain A = ( Jj^> r )5' 1 ] ■ ■ 1 j Jtj^' r )S 
Moreover, the number of shares to invest into asset i is given by ~- {*' r } = J^-. Thus A coincides 

' & J S % (t,r) OSi 

with the classical 'delta hedge'. 

Example 4.8. As in Example 11.1 1 suppose that R is the moving average cHDD process modelled 
as a geometric Brownian motion, and assume that there exists one tradable correlated risky asset. 
More precisely let d = 2, k = m = 1, p = ( a 2 ), (3 = ( ft ft ) with a 2 ,ft,ft G R\ {0}. 
Then 

A /. \ 9p(t,r) ft 

A(t ' r) = " a2 Or # + 

Example 4.9. Applying our results to Example 11.21 we have to take m = 2, k = 2 and d = 3. 
Hence 

.° .° )' P=(b b o)' 

V 72 73 74 / V ft ft u y 
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With a simple minimum square calculation we compute 

n C{ t, r )[VrP(t,r)p(t,r)] = ( li^p(t,r) + 72 ^p(t,r) Js^P&r) ) 
Equation (|l(jp applied to our example produces the following Delta hedge for (t, r) £ [0, T] x R 2 

A(*,r) = ( -^K*,r) + (|f - *)&p(t,r) -£&p(t,r) ) 

where ri represents the crude oil and r 2 the kerosene variable. If 74 = then we have a perfect 
hedge and if 73 = 0, then the price of heating oil doesn't play a role in the hedge, as one would 
expect. 



5 Pricing by marginal utility 

Suppose there is no exchange and the derivative F(Rt) is sold over-the-counter. What is a 
reasonable price a seller could ask for the derivative? The indifference price seems to be a 
natural candidate, though it has the disadvantage that the price of a single derivative depends 
on the total quantity sold, i.e. the indifference price is non-linear. For instance the indifference 
price of 2 x F(Rt) does not equal twice the indifference price of F(Rt). In order to obtain a 
linear version one may take the limit of the indifference price as the quantity converges to 0. 
The object thus derived is the indifference price for a vanishing amount of derivatives, and it 
is therefore called marginal utility price (MUP). Having to pay the MUP for each derivative an 
investor is indifferent between buying and not buying an infinitesimal amount of the derivative. 

We continue requiring (R1)-(R3) to be satisfied. We update the notation and, for q E K and 
(t, r) 6 [0, T] x M m define by p(t, r, q) the indifference price of q units of F(R t rf), i.e. p(t, r, q) is 
the unique real satisfying 

sup{EU{v + G^' r + qF(R^f ) - p(t, r, q)} = sup{EU{v + G% t,r )}. 
A A 

The price of one unit is equal to SiLl^il^ (g ^ 0), and the MUP is defined by 

MUP(t,r) = -^p{t,r,q)\ q=(j . 
Recall that p(t, r, q) = Y t l ' r - y/' r ' 9 , where (y*» r >«, Z 1 ™) is the solution of the BSDE 
Y*™ = qF(R t r f) - f Z^dWu - [ T f(u, R$f, &™)du, s E [t, T]. 

J s J s 

Naming ^(q) = qF(Rj!'), then clearly £(q) is a globally bounded differentiable Lipschitz 
function (with bounded derivatives). The boundedness of £ is trivial since F is bounded and we 
are only interested in the differentiability of the process with relation to q in a neighborhood 
of zero. And so, due to the boundedness of F and the quadratic growth hypothesis for / 
the conditions of Theorem 16.81 are satisfied. Hence, the process Y t,r ' q is continuous in t and 
continuously differentiable in q. 

Writing the BSDE differentiated with respect to q gives 

|-F S '^ = F(Rf) - £ ^Z^dW u - £ V 2 /(n, B*f, Z^^-Z^du, s G [t, T]. 
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Setting q = and renaming the processes for ease of notation we obtain 

U 1 / = F(R t rf) - f V s dW s - C V z f(s, Rl' r , Zf) ■ V s ds. (11) 

J s J s 

As an end product of these calculations we obtain the following explicit formula for the (MUP) 
of our derivative. 

Theorem 5.1. The explicit formula for the Marginal Utility Price of the derivative F(Rt) is 
given by 

MUP(t,r) = U*' r , 
where U t ' r is the first component of the solution pair of the BSDE 

= F(R*f) - f V s dW s - f V z f(s, RY, Z*/) ■ V s ds. (12) 

J s J s 

Remark 5.2. Note that by performing a Girsanov change of probability measure to the one 
making the process W = W + f Q V z f(s, RY , Zt' r )ds a Brownian motion, solving (fT2|) reduces 
to taking conditional expectations with respect to the underlying filtration. Hence, denoting by 
£(■) the stochastic exponential operator, we can represent the marginal utility price explicitly 
by the following expression 

MUP(t,r) = E^^j^VJis^l'^ZlndWs^FiRtf) 

6 Some mathematical tools: smoothness of quadratic FBSDE 

6.1 Moment estimates for BSDE with random Lipschitz condition 

In the following we provide moment estimates for BSDE with generators that satisfy Lipschitz 
conditions with random bounds for the slopes. More precisely, we assume that for our generator 
/ : $7 x [0, T] x IR^ — > R there exists an ]R + -valued predictable process H such that for all 
(u, t,z) EQx [0, T] x R d we have 

\f(u,t,z)-f(u,t,z')\<H t \z-z'\. (13) 

We will assume that H is such that the stochastic integral j HdB with respect to a Brownian 
motion B is a so-called BMO martingale. Recall that J HdB is a BMO martingale (we also say 
it belongs to BMO) if and only if there exists a constant C G M+ independent of u> such that 
for all stopping times r with values in [0, T\ we have 



E 



T 



'Ft 



< C, a.s. (14) 



We refer to |Kaz94j for basic information about BMO martingales. We will abuse the definition 
and refer to the smallest C 6 M+ that satisfies inequality (|14|) as the BMO norm of H. 
Throughout let W be a cf-dimensional Brownian motion. Consider the BSDE 

Y t = Z- [ Z s dW s + / f(s, Z s )ds, < t < T, (15) 
Jt Jt 

where £ is a bounded ^r-measurable random variable, and / satisfies (|13p relative to a pre- 
dictable H with finite BMO norm. 
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We refer to [BC07] for sufficient criteria for the existence of solutions of such BSDEs. 

The moment estimate we shall give next will be needed later for establishing smoothness 
of the solution of the quadratic BSDE with respect to the parameters the terminal condition 
depends on. 

Lemma 6.1. Suppose that for all (3>\ we have \f(s,0)\ds G L^{P). Let p > 1. Then there 
exist constants q > 1 and C > 0, depending only on p, T , and the BMO-norm of H , such that 
we have 

i 















E 


sup \Y t \ 2 P 


+ E 




< c (e 


~\Z\ 2pq + ([ T \f(s,0)\ds) 2pi1 ' 




L te[o,T] 








1 Jo J 



Proof. This follows from Corollary 3.4 in |BC07j by a straightforward generalization to the 
multidimensional case considered here. It can also be shown with the method used in the proof 
of Theorem 5.1 in |AIDR07| . □ 

6.2 Differentiability of quadratic FBSDE 

Consider now a FBSDE of the form 

X x s = x + j*b(s,X*)d S + ftp(s,X*)dW s , 

Yf = F(X$)- jf ZfdW, + jf f(s,X5,Zf)ds, 

where b : [0, T) x R m — > R m and p : [0, T] x R m -> R mxd and W is the d- dimensional Brownian 
motion of the preceding subsection. Note that p is a n x d matrix. We will denote its transpose 
by p*. The generator of the backward part is assumed to be a 7 3 (.Fj)®;B(M' m )(g>,8(]R rf )-measurable 
process / : f2 x [0, T] x R m x R d — > R such that there exists a constant M & R + such that for 
all (t, x, z) G [0, T) x R m x R d we have 

\f(t,x,z)\ < M(l + \z\ 2 ) a.s. (17) 

Here V{Tt) denotes the u-field of predictable sets with respect to the filtration [T-t). Moreover 
we assume that 

/ is differentiable in x and z and , . 

\V z f{t,x,z)\ < M(l + \z\) for all (i, x, z) G [0, T] x R m x R d a.s. ^ ' 



We will give sufficient conditions for the process Y x in the solution of the FBSDE (]16p to be 
differentiable in x. A further assumption we need is that the coefficients of the forward equation 
belong to the function space B mxd and B mxl respectively (see Definition 14. 1 1) . To simplify 
notation, to the pair (b, p) of coefficient functions we associate the second order differential 
operator £ = ZZi fc(0& + \ E™=i[PP%Q^- 

We will assume that the coefficients of the forward equation (|16p satisfy 

(Dl) p G B mxd , b G B mxl , 

and that 

(D2) F : R m -> R is a twice differentiable function such that VF • p G B lxd and £F G B lxl . 

It is known that the conditions (Dl) and (D2) ensure that X x is differentiable in x and the 
difference quotients can be nicely controlled. For the convenience of the reader we quote a 
standard result which will be needed later. Denote by the unit vector in R m in the direction 
of coordinate i, 1 < i < m. 
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Lemma 6.2. Suppose (Dl) and (D2) are satisfied. For all x G M. m , h 7^ and i G {1, . . . ,m}, 
let C xXi = \(F(X^ +hei ) - F(X%)). Then for every p > 1 there exists a C > 0, dependent only 
on p and the bounds of b, p, F and its derivatives, such that for all x, x' G K m and h, h! 7^ 0, 



E 



IC 



x.h.i 



c 



x',h',i\2p 



< c(\ x -x'\ 2 + \h-h'\y. 



(19) 



Proof. Note that by Ito's formula F{Xf) = F(X§) + jj VF(X*) • p(s, X* )dW s + Jj| £Fds. Thus 
F(Xf) is a diffusion with coefficients p(s, rr) = VF(x) ■ p(s, x) and b(s, x) = YliLi h(s, ■ 



, 8F(x) 
dxi 



+ 



E=lft)'(*-2)3, M € M x Km - % (D2) we have p G B lxd an d b G B lxl . 
Therefore, by using standard results on stochastic flows (see Lemma 4.6.3 in |Kun9()| ), we 
obtain the result. □ 



Notice that since F is bounded and growth condition (|TTj) holds, there exists a unique solution 
(Y X ,Z X ) G <S°°(R) <g> H 2 (R d ) of the BSDE in (JTgJ) for all x G M d . One can even show that we 
may choose the family (Y x ) x ^m such that it is continuous in x. 

Lemma 6.3. Let (Dl), (D2), \1T\ ) and A18\) be satisfied, and assume that F is bounded and that 
there exists a constant K G M+ such that for all t G [0, T], x, x' G and z G R d 



\f(t,x,z)-f(t,x',z)\<K(l + \z\)\x-x'\. 
Then for all p > 1 there exists a constant C G M+ such that for all x, x' G 

,'|2p 



E sup \Y X -Y x '\ 2 p < C 
te[o,T] 



x — x 



E 



T 



Z x '\ 2 dt 







< C\x — x 



i\2p 



(20) 

(21) 
(22) 



In particular, Kolmogorov's continuity criterion implies that there exists a measurable process 
Y : U x [0/T] x R m such that (t, x) i-> is continuous for a. a. u; and for all (t, x) G [0, T] x R m 
a.s. 



we have Y x 



Proof. For a G M, let x(y) = e ° y - By applying Ito's formula to x(Y x ) an d using standard 
arguments one can show that f Q Z x dW G BMO with the BMO norm depending only on the 
bound of F and the growth constant of / in z. 

For all x,x' G R m let U t = Y t x - Y t x ' , V t = Zf - Zf' and C = F{X X ) - F{X X '). We use a 
line integral transformation in order to show that U x can be seen as a BSDE with generator 
satisfying a Lipschitz condition of the type (fT3j) . Define Jt = V x f{t,X® — $(X£ — Xf), Z^dfi 
and H t = Jq 1 V z f(t, Xf, Z x ' - 0(Z? - Z?'))d$ and observe that 



fv s dW s + f T (f(s,X x ,Z x ) - f(s,X x ',Z x )) + (f(s,X x ',Z x ) - f(s,X x ',Z x '))ds 
Jt Jt 

[ T V s dW s + f T {J s {X x - Xf) + H s V s )ds. 
Jt Jt 



The moment estimate of Lemma 16 . 1 1 applied to the pair (U, V) leads to 











E 


' sup \U t \ 2p 


+ E 


(f iv.i'*)' 




L te[o,T] 







< C E 



ICI 



2pq 



+ 



T 



\J S {X X -X x ')\ds) 2 ™ 
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for some constants C > and q > 1. By (|20p we have V x f{t,x,z) < K(l + |z|), and hence 

£(^ T |J S (X* - X*')|ds) 2p(? < K 2pq + l^'l) 2 ^) 2pg ) 2 " ^'I 2 ^) 2w ) 2 • 

Lemma 16. II implies that E(J^ (1 + |Z x |) 2 (is) 2p9 is bounded, and by standard results on moment 
estimates of SDEs we have E(Jq \Xf - Xf\ 2 ds) 2pq < C'\x - x'\ ipq for some constant C G M+ 
(see Theorem 3.2 in |Kun04j ). Moreover, the Lipschitz property of F guarantees that there 
exists a constant C" such that E\(\ 2pq < C"\x — x'\ 2pq , and hence the desired result follows. □ 

The following theorem guarantees pathwise continuous differentiability of an appropriate 
modification of the solution process. 

Theorem 6.4. Let (Dl), (D2), \17\ ) and $18\) be satisfied, and suppose that F is bounded and f 
satisfies 120\) . Besides suppose that V z f is globally Lipschitz continuous in (x,z) and that V Xi f 
satisfies for all t £ [0, T], x, x' € and z, z' G M. d 

\V x f(t,x,z) - V x f(t,x',z')\ < K(l + \z\ + \z'\)(\x - x'\ + \z- z'\). (23) 

Then there exists a function n x [0,T] x R m -» M. m+1+d , (w,t,x) ■-> (Xf,Y t x ,Z?)(w), such that 
for almost all uj, Xf and Y x are continuous in t and continuously differ entiable in x, and for 
all x, {Xf ,Y t x , Zf) is a solution of FBSDE [TS\) . Moreover, there exists a process V X Z X € 7i 2 
such that the pair (V X Y X ,V X Z X ) solves the BSDE 

V X Y X = V x F(X x )V x X x -f t T V x Z x dW s 

+ f[V x f(s,X x ,Z x )V x X x + VJ(s,X x ,Z x )V x Z x }ds. 



We will use Kolmogorov's Lemma in order to prove the theorem. Let x € M. m . For all 
h + 0, let A x t > h = i(Xf+^ - Xf), Uf> h = \{Y x ^ h - Y x ), V? h = \{Zf +h ^ - Zf) and 
( x,h = r(£(x + hei) — £(x)). We need the following estimates. 

Lemma 6.5. For all p > 1, x, x' £ M m , h,h' ^ we have with some constant C 



E 



sup \U t — U t 

te[o,T] 



<C(\x-x'\ 2 + \h-ti\ 2 ) p . (25) 



Proof Let p > 1. Note that for all h ^ 

U?' h = C' h ~ f Vt h dW s + J T ~[f(s,Xf +he *,Zf+ h et) - f(s,Xf,Zf)]ds. 

We use a line integral transformation in order to show that U x,h can be seen as a BSDE with 
random Lipschitz bound. To this end define two (^t)-adapted processes by 

A x t > h = [\ x f{t,X? + 0(Xf +h «-X?),Zf)d0, 
Jo 

lf' h = [ V z f{t,Xf +he %Zf + $(Zf +he * - Zf))dd. 
Jo 



Then 



i[/(t, Xf +he \ Zf +he *) - /(*, Xf, Zf)} = A x t > h A x t > h + lf> h V?> h . 
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The growth condition jTS]) implies that \I x > h \ < M(l + \Z X \ + \Z x+he *\), and hence f Q I x > h dW G 
BMO. Thus we obtain a BSDE with generator satisfying condition (|13j) . 

Now let G M m and h,h' + 0. Then the difference (U x ' h - U x '> h ',V x > h - V x '' h ') solves 
again a BSDE with generator of the type (|13fl . namely 



fx' \h' 



z, dW, 



pT 

- J (I x ' h z s + (I x > h - ifV)V x ' h + Af h Af h - A x '> h 'A x '' h ')ds. 
Therefore Lemma 16 . 1 1 yields for q > 1 



E 



\r T x,h TT x'h! \2n 

sup \U t ' - U t I p 
*e[o,T] 



< C\E 



\^£X,h ^x' ,h' |2pg 



+E 



{\A x ' h A x / - A x ^ h ' A x ^ h '\ + \I\ 



rx,h 



To treat the first term on the right hand side, use LemmaEJto see that E[\C x ' h -( x '' h '\ 2pq ] « < 
C(|x-x / | 2 + |/i-/i'| 2 )f. 

For the second term, recall that V z / is Lipschitz continuous, say with Lipschitz constant 
L G R + . We therefore have for any t G [0, T] 



I jx,h jx' ,h' I 



< L(\X x ' h - X x > h | + \Zf - Zf | + | Z* +ftei - Zf +,r 61 1 ) . 



Now Cauchy-Schwarz' inequality leads to 

rT 



E 



< E 



jX,h_jX,h ^yx ,h 







\I x > h - I x '> h '\ 2 ds 



2pq 



2pq ^ i 
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So Lemma 16.31 and Lemma 4.5.6 in [Kun90j imply with some constant C 



E 



2pq 



if> h -if> h '\ 2 ds) 2q < c(\ x - xf + \h - hfy . 



The term E 



C\V x ' h \ 2 ds 



2pq 



is seen to be bounded by an appeal to Lemma [67 



J^\A x ' h A x ' h -A x '' h 'Af' h '\dsY" q <C(\x-x'\ 2 + \h-h'\ 2 )P 



2pq 



It remains to show that E 

For this we separately estimate the two summands on the right hand side of the following 
inequality 

| ^x,h /\x,h A x>,h' A x',h'\ <^ | ^x,h | >£x,h ^x',h' I _|_ j ^x' ,h' | 1 ^x,h ^a;',/i'| 

First note that due to (f!ZU|) we have for some constants C\ , C2 . . . 

£ \Af h \\Af h - A x '> h '\ds < d Q^(l + \Z x \) 2 ds\ 2 ^£ \Af h - A x '' h '\ 2 ds\ * , 
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which implies, together with Lemma 16. II and standard estimates of differences of the A x,h (see 
Theorem 3.3 in |Kun04] ). 




2p<?\ 2 



< C 3 (\x-x'\ 2 + \h-h'\ z ) pq . 



Secondly, from (jz3|) we obtain 



/ \A x '> h '\\Af h - A x s '> h '\ds 
Jo 

< C 4 / (1 + |Z*| + |zf |) (p£? - + |X* +fte * - X? +h,ei \ + |Zf - Z x '\) ds 



o 



and hence, with Lemma 16.31 and moment estimates for X 1 

iX ' h -A x s '' h '\ds 




[ J [\Xf - Xf | + \X x+he * - X x ' +h ' e * | + \Z X - Zf | 

< C 6 (\x-x'\ 2 + \h-h'\ 2 ) pq . 
Combining the estimates just derived, we conclude 



2p<A 2 



E 



x,h ttX ,h' \2p 



sup \U t — U t 

te[o,T] 



< C{\x-x'\ 2 + \h-h'\ 2 Y 



This completes the proof of the lemma. □ 

Proof of Theorem \6.4\ Note that by Lemma 16.31 we may assume that (t,x) i— ► Y x is continuous 
for all lo. Then U x,h has continuous paths for all x G W 11 and h ^ 0. 

Let Q be the collection of all pairs (x, /i) where x is a vector of dyadic rationals in M m and 
ft / a dyadic rational in R. Since inequality (|25p is valid, Kolmogorov's lemma implies that 
there exists a null set N such that for all to G A^ c the function Q 3 (x, /i) i— > U x,h can be uniquely 
extended to a continuous function from R m+1 into the space of continuous functions endowed 
with the sup norm (see Thm 73, Ch. IV, |Pro04j ). Such a null set N can be chosen for any 
direction i in which we differentiate, and hence there exists a modification of Y x such that for 
all t the mapping x i— > Y/^ possesses continuous partial derivatives. 

Finally it is straightforward to show that the derivative V X Y X together with a process V X Z X , 
defined as an U 2 limit of the processes V x,h as h — > 0, solve the BSDE (|24p . □ 



6.3 The Markov property of FBSDE 

The forward part of our FBSDE (|16p is solved by a time inhomogeneous Markov process. We 
will now investigate the consequences of this fact in more detail. Let us fix an initial time 
t £ [0, T), as well as an initial state x to be taken by our forward process at this time. Then, 
conditioned on taking the value x at time t, the forward process satisfies the SDE 

Xl> x = x+ f b(r, X l r < x )dr + r p(r,X}> x )dWr, (26) 
Jt Jt 
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where x G R m and s G [t,T]. We will assume that the coefficients satisfy a growth and a 
Lipschitz condition. More precisely, assume that there exists a constant C G R+ such that for 
all x, x' G R m and t G [0, T] 

\b(t,x)-b(t,x')\ + \p(t,x)-p(t,x')\ < C(\x-x'\), 

\b(t,x)\ + \p(t,x)\ < C(l + \x\). W 

Condition (|2Tj) guarantees that there exists a unique solution of (|26p . It moreover implies that 
is Malliavin differentiable and that its Malliavin gradient has a representation involving, 
for (t, x) fixed, the global flow on the space of nonsingular linear operators & ,x on R m defined 
by the equation 

&f = 1r«» + J' V x b(u,X^)&>*du + J" V x p(u,X^)^dW u , s > t. 

Here V x 6 and V x p describe the gradients of b resp. p existing in the weak sense under (|27|) . 
iRm the m x m unit matrix. The Malliavin gradient is then given by the formula (see Nualart 
|Mua95| . p. 126) 

DoX** = $*/($f r 1 X 1 /), t<$<s. (28) 

With the Markov process X t,x starting at time t in x in mind, we now consider BSDE of the 
form 

Y s t ^ x = F(X t /)- Z t r > x dW r + f^X^^Z^dr. (29) 

J s J s 

In accordance with Section [21 we now assume that the generator is a deterministic Borel mea- 
surable function / : [0, T] x R m x M. — > M. Again we assume that / is differentiable in (x, z) 
and that there exists a constant M G M+ such that for all (t, x, z) G [0, T] x R m x M. d we have 

\f(t,x,z)\ < M(l + \z\ 2 ) a.s. and \V z f{t,x,z)\<M(l + \z\) a.s. (30) 

for all (t, z) G [0, T] x R m . If F is bounded, then it follows from Theorem 2.3 and 2.6 in [KobOO] 
that there exists a unique solution (Y^, Z l ' x ) G S°°(IR) ® H 2 {R d ) of the BSDE The next 

result states that the solution of the BSDE is already determined by the forward process X l ' x . 
In order to formulate it, for all m G N we denote by T> rn the u-algebra on R m generated by the 



family of functions M. m B x i— ► E fT ip(s, Xt' x )ds, where t G [0, T] and (p : [0, T] 



is 



bounded and continuous. 



Theorem 6.6. Let F : R m —> R 6e a bounded Borel function, suppose that f satisfies 130\) and 
the coefficients of the forward diffusion \21\j . Suppose that there exist functions f n : [0,T] xR m x 
R d — > R, globally Lipschitz continuous in (x, z), such that for almost all oo and for all compact 
sets K C R m x R d the sequence f n converges to f uniformly on [0, T] x K. Then there exist two 
B[0, T] <g) D m - and B[0, T] <g> T> m -measurable deterministic functions u and v on [0, T] x R m such 
that 

Y^ x = u(s, Xtn and Z^ x = v(s, X^s, X^), (31) 
for P <g) A-a.a. (w,s) efix [t,T]. 

Proof. Let / n be Lipschitz continuous in (x, z) such that / n converges locally uniformly on 
R+ x R m x R d . Let (t, x) G [0, T] x R m and denote by (Y n , Z n ) = {{Y n ) l ' x , (Z n ) l ' x ) the solution 
of the BSDE with generator f n and terminal condition £ = F(Xj, x ). It follows from Theorem 
2.8 in [KobOOj that Y n converges to Y f > x in H°°(R), and Z n converges to Z*' x in H 2 (R d ). By 
taking a subsequence if necessary, we may assume that Z n converges to Z t,x a.s. on Q x [0, T]. 
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According to Theorem 4.1 in |KKPQ97| , there exist B[0, T]®V m - and B[0, T]® Immeasurable 
deterministic functions u n (t,x) and v n (t,x) that satisfy the representations Y" s n = u n (s,Xl' x ) 
and Zg = v n (s, X t s ' x )p(s, Xt' x ) for all s S [t,T] a.s. Now define 

u(t, x) = liminf u n (t, x) and v(t, x) = liminf v n (t, x). 

n n 

Then Yt x = u{s, Xt' x ) and Z l f = v{s, X t / X )p{s, X l / X ), a.s. □ 

By combining Theorem 16.61 with Theorem 16.41 we obtain a representation of the control 
process Z t,x in terms of the derivative of Y t,x with respect to x. 

Theorem 6.7. Suppose that the assumptions of Theorem \6.(ft are satisfied. Besides assume 
that V z / is globally Lipschitz continuous, that $2Q\) and 123\) are satisfied, and further that the 
forward coefficients satisfy the stronger conditions (Dl) and (D2). Then u(t,x) is differentiable 
in x for a. a. t £ [0,T]. Moreover, 

Z^ = V x u(t,Xl' x )p(s,Xt' X ), (32) 

for P ® X-a.a. (u,s) € Q x [t,T]. 

Proof. Recall that Xs' x is Malliavin differentiable and that the assumptions of Lemma 16.31 are 
satisfied. Equation (f2T|) implies that x i— > u(t, x) = Y^' x is Lipschitz continuous. Therefore 
Yg' x = u(s,X s ,x ) is Malliavin differentiable (see Proposition 1.2.2 |Nua95j ) . By Theorem 16.41 
u(t, x) is differentiable in x, and by the chain rule we have D$Ys' x = V x u(s, Xs X )D$xt' x . Since 
due to (|28p D s Xt' x = p(s,Xt' x ) and Zt' x = D s Yg ,x (the later following f.ex. from Lemma 5.1 in 
|EKPQ97| ), Theorem [62] implies ([32]). □ 



6.4 Differentiability of Quadratic BSDE with parameterized terminal condi- 
tion 

For this subsection we pass to a more abstract parameter dependence of the solution of a BSDE 
than studied above in a pair of forward and backward SDE. We consider the BSDE 

Yt = - [ T Z x dW s + f T f(s, Z x )ds, tE[0,T],xe R m . (33) 

Jt Jt 

Throughout we assume that 

(El) 1R" 1 9ih £ f is a bounded random field which as a function of x is differentiable 
with bounded partial derivatives; V£(x) is also Lipschitz in x; also f(t,0) is (Pt)— adapted 
and satisfies fit, 0) € LP for all p > 1. 

(E2) there exists M £ M+ such that \f(t, z)\ < M(l + |^| 2 ) a.s.; / is differentiable in z such that 
\V z f{t,z)\ < M(l + \z\) for all (t,z) G [0, T] x R d a.s. 

(E3) for all x G R m , h ^ and i E {1, . . . ,m}, let Q x,h,i = + hei) - £(x)). Then for every 

p > 1 there exists a C > 0, dependent only on p, such that for all x, x' £ M. m and h, h! ^ 0, 

<C(\x-x'\ 2 + \h-h'\ 2 ) p . (34) 

Although the terminal condition does not depend on a forward diffusion (see Lemma 16.21 for 
a derivation of (|34|) in a FBSDE setting), Hypothesis (E1)-(E3) allow to apply the methods we 
used in Subsection 16.21 and obtain the following theorem. 



E 



IC 



x.h.'i 



,h',i\2p 
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Theorem 6.8. Let (El), (E2) and (E3) be satisfied. Then there exists a function Q x [0, T] x 
R m — ► R 1+rf , (oj,t,x) i — > (Y x , Z x )(uj), such that for almost all uj , the process Y x is continuous 
in t and continuously differentiable in x, and for all x, (Y X ,Z X ) is a solution of BSDE h33\). 
Moreover, there exists a process V X Z X G 7i 2 (M. mxd ) such that the pair (V X Y X , V X Z X ) solves the 
BSDE 

V x Y t x = V^(x) - / V x Z x dW s + / [V z f(s, Z X )V X Z X ] ds. 
Jt Jt 

Proof. Conditions (El) and (E3) guarantee that the solutions of the BSDE (|33|) exist and 
(Y X ,Z X ) e 5°°(R) ® H 2 (R d ). 

Condition (El), (E2), (E3) and the BMO property of the martingale J Z x dW allow us to 
prove moment estimates that correspond to Lemma I6.1( Lemma 16.31 and Lemma 16.51 Hence a 
simple adaptation of the proof of Theorem 16.41 provides the proof of Theorem 16.81 □ 
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